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Abstract. 

The Dzyaloshinskii-Moriya interaction is shown to result in a canting of spins 
in a single molecule transistor. We predict non-linear transport signatures of this 
effect induced by spin-orbit coupling for the generic case of a molecular dimer. The 
conductance is calculated using a master equation and is found to exhibit a non-trivial 
dependence on the magnitude and direction of an external magnetic field. We show 
how three-terminal transport measurements allow for a determination of the coupling- 
vector characterizing the Dzyaloshinskii-Moriya interaction. In particular, we show 
how its orientation, defining the intramolecular spin chirality, can be probed with 
ferromagnetic electrodes. 

PACS numbers: 73.63.Kv , 85.65.+h, 85.35.-p, 71.70.Ej, 85.75.-d 

Submitted to: Nanotechnology 
1. Introduction 

Antisymmetric exchange interaction first identified by Dzyaloshinskii based on 
symmetry arguments p], and derived by Moriya [2], is fundamental to the understanding 
of weak ferromagnetism in materials like a — Fe2C>3 [3]. Nowadays such spin-orbit 
effects can be addressed in nanoscale experiments where transport plays an important 
role, for instance, for STM on atomic chains or even single magnetic atoms on a 
surface [3] . In particular, the Dzyaloshinskii-Moriya (DM) interaction favours spiralling 
spin-structures in systems with broken inversion symmetry, which is generically the case 
on a surface [5]. Intrinsic effects of spin-orbit coupling in magnetic molecules have also 
been investigated in three terminal devices [HI [7J [8], focusing on easy-axis magnetic 
anisotropy and quantum-spin tunnelling effects [9j [10] . Whereas the latter are induced 
in second order in the spin-orbit (SO) interaction, the DM interaction already arises 
in linear order. In particular, DM interaction is known to have a large effect on the 
magnetism of dimer molecules [11] and plays an important role in the single molecular 
magnets Mn 6 p], Mn 12 M, EE] and Vi 5 [H EH] • 

For the electronic spectroscopy of single molecule transistors (SMTs) [17] the spin- 
orbit interaction is of special importance, since it may generate a specific electric or 
magnetic field dependence of the transport [18]. This facilitates the experimental 
discrimination between the molecule's contribution to the electric current and that of 
possible spurious particles in a junction. It also leads to a violation of spin-selection 
rules, which clearly affects the non-equilibrium occupations of the magnetic molecular 
states [9]. Since the Dzyaloshinskii-Moriya interaction relies fundamentally on broken 
spatial inversion symmetry, which is typically present in a molecular junction, it is of 
particular interest to identify its characteristic transport signatures. 

Dzyaloshinskii-Moriya interaction arises in its most elementary form in a molecular 
dimer [11] where two localized electrons interact via Coulomb and hybridization terms, 
in combination with local spin-orbit scattering into excited orbitals. As a result, an 
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effective spin-spin interaction arises between the spin 1/2 localized at the monomers 
labelled by i — 1, 2 (see figure [2]): 

H ea = JSi ■ S 2 + D ■ (S 1 x S 2 ) + B ■ (Si + S 2 ) (1) 

The isotropic exchange typically dominates, resulting in a singlet and triplet split by 
J. In three terminal molecular devices [TH [20] such a splitting in the meV range has 
been reported. The antisymmetric exchange interaction is specified by the DM vector 
D which is an intrinsic property of the molecule. This interaction typically presents a 
weak perturbation which mixes the two multiplets of different total spin length, S — 
and 1 here, resulting in small splittings ~ |D| 2 / J <C J. The magnetic field B (in energy 
units, where gfiB = 1, assuming an isotropic g = 2) enhances this mixing by bringing 
one triplet component close to the singlet. In general, this results in an energy level 
anticrossing at |B| « J (see figure 1(a)). Interestingly, this mixing strongly depends on 
the orientation of the field B relative to the DM vector D. To see this, we quantize the 
spin along the magnetic field, i.e. B is chosen along the ^-direction. For B || D the singlet 
state does not mix with the lowest triplet and the states can cross, since in this case 
the Zeeman and DM terms commute (they contain three different spin-components). 
In contrast, for B _L D the singlet and S x — ±1 triplet states are mixed. As a result, 
when sweeping the magnetic field there is a smooth adiabatic change of the ground 
state magnetization instead of an instantaneous jump at the crossing point |B| « J (see 
figure 1(b)). Furthermore, the relative orientation of the two spins shows an interesting 
field dependence: For example, in case of the ground state the antiparallel spins are 
slightly canted (left of figure 1(b)) at zero field and rotate towards each other with 
increasing field, until they are almost perpendicular at the anticrossing point B « J. 
For higher fields B > J they become aligned by the external field (right of figure 1(b) ). 




Figure 1. (a) Eigenenergies and (b) spin projection onto the magnetic field B, 
which points along the z-axis and is perpendicular to the DM vector D = (0.2 J, 0,0). 
The red (green) curves and arrows in (a) and (b) refer to the ground (first excited) 
state. The arrows in (b) indicate the relative orientation of the monomer spins. The 
isotropic exchange J is antiferromagnetic, resulting in near antiparallel orientation of 
the monomer spins in the ground state at low fields (red curves). 
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Figure 2. Dimer with localized spins, which are slightly canted due to interplay 
between DM interaction and isotropic exchange coupling J. The dimer is probed in 
a three terminal junction (Vb- Bias voltage, Vg- Gate voltage). Tunnelling is only 
possible between monomers and from one monomer to the adjacent lead. The grey 
vertical lines indicate the spin orientation without DM interaction. 

The expectation value (Si x S2) quantifies the "canting" of the spins. In particular, its 
sign equals the chirality characterizing the intra-molecular spin structure. 

In this work we predict the single electron transport signatures of the most 
elementary realization of DM interaction in a molecular dimer placed in a three terminal 
junction (see figure [2]). In section [2] we introduce a microscopic model which captures the 
basic physics discussed above and comment on the calculation of the transport current 
using a master equation. In section 13.11 we present the transport current as function 
of both bias voltage and magnetic field. We show that the DM interaction leads to 
characteristic transport features which depend strongly on the gate voltage. We propose 
magnetic field measurements which allow for a determination of the magnitude of the 
D- vector and the axis parallel to it. Finally in section 13.21 we present the non-linear 
conductance for polarized electrodes with antiparallel magnetizations. We show that 
due to the DM interaction the transport is sensitive to the reversal of both electrodes' 
polarizations. We explain how this additional measurement allows the polarity of the 
D- vector to be determined as well. 

2. Model 

To predict the basic signatures of transport through a "molecular double-dot" as 
sketched in figure [2] including the effect of DM interaction a customary spin-Hamiltonian 
is insufficient: multiple charged states as well as orbital degrees of freedom of the 
monomers need to be accounted for. Instead one has to consider the underlying many- 
body molecular Hamiltonian [2] from which the Dzyaloshinskii-Moriya interaction ([1]) 
derives: H m = H Q + Hi + H so + + Hz- In particular, both ground {i = 1,2) 
as well as excited orbitals of the monomers (i = 1', 2') have to be included through 
H Q = J2ia e i n i° with energy e^. The spin-orbit scattering with amplitude A« between 
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these orbital states on monomer % = 1,2 reads: if so = ^\ . a a i\d\ a hj • Saa'dj a >- Here 
dicj and n.j CT = d\ a di a denote the electron and number operator of orbital i with spin 
a and s aa t are the spin | matrices. We assume for simplicity no orbital degeneracy 
implying that the orbital angular momentum is quenched on each monomer (1& = 0). 
The off-diagonal matrix elements of the orbital angular momentum 1, which parametrize 
the local molecular orbital structure, are then imaginary and have a crucial effect on 
the spin-orientation through the SO interaction. The magnetic field tunes this effect 
through the Zeeman term H z = B • (|l + s) , which acts on the spin s = ^ a a , d\ a s aa / d ia > 
and the orbital momentum 1 = £\ . a d ia lijdi a . Characteristic for a covalently bonded 
dimer is the hybridization = ■ a d\ a tijdj a + h.c\ between the monomers and 
it couples their excitations [2TJ [22]. The strong local electron-electron interaction, 
H i = \Y.o,o>Y.ijki v ijki d \ (T d) (T ,d la ,d ka , plays a crucial role in several ways. For the 
transport, local repulsion results in Coulomb blockade of tunnelling and the transport is 
dominated by single electron tunnelling (see below). Furthermore, the direct (Hund's) 
exchange {v^i = Jh), combined with antiferromagnetic superexchange induced by the 
hybridization and the local repulsion, gives rise to the isotropic exchange interaction 
of the spins in (0Q) with J = Jh + jj and U = van — Vijij. The Hamiltonian H m is 
diagonalized for different values of the electron charge N of the dimer (N = 0, 1, 2). 

The above model incorporates the processes which generate the Dzyaloshinskii- 
Moriya interaction in a way analogous to isotropic exchange. In linear order in spin- 
orbit coupling one can derive from this the effective spin Hamiltonian (CQ) for N = 2 
with 

D = 2 — h uiaai' -(1^2) (2) 



e v - ei V U 

where 1^2 denotes the same term with the roles of the monomers interchanged. 
Here virtual processes involving tunnelling and spin-orbit scattering gives rise to the 
first contribution (see figure 3(a)). In addition, Coulomb excitation of one monomer 
by the other through %j-fc, sketched in figure 3(b), induces an interaction of the 
same form, contributing a second term. Equation §2§ explicitly shows that a non- 
vanishing Dzyaloshinskii-Moriya vector is a joint effect of the asymmetries in the 
monomer excitation spectra, interaction matrix elements and spin-orbit couplings. For 
definiteness we introduce the required asymmetry through the spin-orbit coupling on 
the monomers, Ai = 2A2 = 8meV, keeping other parameters symmetric. Chemically, 
this asymmetry may be controlled by an appropriate choice of coordination ligands. For 
molecular dimers a hierarchy of energy scales with U, > J, A is expected. The results 
discussed here, based on the parameters in table [lj are representative of this order. The 
corresponding energies in the effective spin-Hamiltonian are 

J = 4 meV D = -(0.32 meV, 0, 0) (3) 

i.e. the effective exchange J is antiferromagnetic. By the choice of parameters processes 
indicated in figure 3(b) do not contribute. We checked that they lead to the qualitatively 
similar results. 
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Table 1. Model parameters. Double orbital occupation or occupation of an excited 
orbital is only virtually possible for these values. The resulting virtual spin scattering 
gives rise to the DM interaction. 

(a) (b) 

One-particle interaction Two-particle interaction 

ej OmeV i 6(1,2) 
e v 50meV i'e(l',2') 
tij 10 meV j ^ i' 

(i,0,0) ie(l,2) 



v Hii 125 meV i 6 (1, 1 , 2, 2 ) 
Uji'jj/ 125 meV i 6 (1,2) 
Vijij 25 meV j ^ i' 



Transport in three terminal molecular junctions [T7] can be modeled using a 
tunnel Hamiltonian approach, parametrized by tunnel rate constants T r = T for 
tunnelling to electrode r = L, R. For simplicity, these are assumed equal for the 
ground and excited orbital. Since the basic transport signatures of Dzyaloshinskii- 
Moriya have, to our knowledge, not been discussed before, we focus on the basic single 
electron tunnelling transport regime (SET). In our calculation we assume an electron 
temperature T = 0.02 meV in order to resolve the energy scale set by D in ([3]). In 
the SET regime the stationary occupations of the molecular many-electron states and 
transport current can be calculated from a master equation |9j. Some comments are 
in order. First, off-diagonal elements of the density matrix in the molecular energy 
eigenbasis can in principle be important here, since neither the spin length S nor its 
projection S z are conserved quantities in tunnelling. These are not accounted for in 
the master equation. However, as we restrict our calculation of the tunnelling rates 
to the first order in the tunnel coupling and the energy spectrum contains no quasi- 
degeneracies, these elements can be neglected. We have also checked this explicitly 
by comparing with results from a generalized master (or kinetic) equation accounting 
for non-diagonal elements. Secondly, we note that in magnetic molecules long spin-life 
times have been reported [23]. For this reason we focus mainly on calculations neglecting 
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Figure 3. Virtual processes leading to Dzyaloshinskii-Moriya interaction. (a) 
superexchange and (b) direct exchange. 
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any magnetic relaxation unrelated to transport and find signatures of non-equilibrium 
magnetic states. Still, to appreciate the possible influence of relaxation process (due e.g. 
coupling to junction phonons etc.), we have for each calculation performed a worst-case 
estimate calculation in which all excited states relax to all lower lying states within the 
same charge sector with a rate up to the tunnel rate (We assume a quadratic energy 
dependence of the density of states for these energy loss processes). We thus identify 
transport signatures which are sensitive to relaxation and those which are robust against 
it. 

3. Transport signatures of DM interaction 

3.1. Normal electrodes 

For orientation, we show in figure H] the expected resonance lines for a dimer without 
DM interaction in a so-called stability diagram [23], i.e. oII/oIVb plotted versus the 
bias, Vb, and gate voltage, Vq- We have explored this map in detail and identified 
significant transport signatures of Dzyaloshinskii-Moriya interaction along several sets 
of transitions, which are marked by I, II, III and IV in figure HI Due to the spin-orbit 
mixing four additional electronic transitions can be observed which violate the spin 




Figure 4. Schematic stability diagram as function of experimentally controllable 
voltages Vb and Vq for low magnetic field B < J and no spin-orbit interaction. 
The grey-shaded area indicate the Coulomb blockade (CB) regime, where no current 
flows. Several SET resonances are marked by the box and the arrows. At these lines 
clear signatures of the Dzyaloshinskii-Moriya interaction develop when the spin-orbit 
interaction is included. 
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Figure 5. Schematic many-body excitation spectra for N = 1 and N = 2 in a 
high magnetic field B z < J relevant to regime I in figure [U The energy scales of the 
hybridization t and isotropic exchange J are indicated. The grey arrows indicate the 
Zeeman splitting of the multiplets. The Dzyaloshinskii-Moriya interaction mixes the 
two low-lying N — 2 states most strongly resulting in energy changes on the order of 
|D| <C J, t. The spin-orbit coupling also leads to spin- forbidden transitions (red dashed 
lines), in addition to the spin-allowed ones (blue solid lines). The shown transitions 
involving the symmetric excited one-electron states are labelled with an asterisk. 



selection rule AS Z — ±|: 

| i) a -> \T +1 ) , | T) a - |T_!> (4) 

where the index a indicates whether the symmetric (regimes I and II) or the 
antisymmetric excited one-electron states (III and IV) are involved. Since the latter 
transitions start from an excited state, they are sensitive to intra-molecular relaxation 
as our calculations bear out. Apart from this the results in all regimes are qualitatively 
very similar. From hereon we focus on the regime I, where the situation is most clear 
since only charge states N = 1 and N = 2 of the dimer are involved. Experimentally, 
this is the most relevant regime since the low voltage Vb ~ J avoids complications due 
to e.g. junction instabilities. In figure [5] we show the relevant transitions between the 
dimer states due to transport. The ground state is labelled with , i.e. indicating 

the intramolecular spin structure. It has mainly singlet character before the anticrossing 
and for high magnetic field closely resembles the pure polarized triplet state. For the 
first excited state |\^), which anticrosses with the ground state, it is the other way 
around (see figure 1(b)). 



We first consider the case where BID. In figure 6(a) we show dl/dVs as the 
bias Vb (along line "A" in figure HJ) and the magnetic field is swept. Without spin- 
orbit coupling, one expects two degenerate transitions at low fields |B| <C J. In 
contrast, two additional weak oII/oIVb peaks (lines [a],[h]) are found with anomalous 
field dependence corresponding to a tripled g-factor since |AS 2 | = § is allowed by the 
spin-orbit interaction. This allows for their correct identification as due to breaking 
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Figure 6. Magneto-conductance maps for regime I: dependence of dl /cIVb and the 
magnetic field magnitude B z for BID. The conductance is normalized to T. (a) 
Overview of all transitions using a logarithmic scale for a voltage sweep along line "A" 
in figure [4] The spin-forbidden transitions are weak and have an anomalous strong 
field dependence (triple g- factor). The isotropic exchange constant J can be extracted 
as indicated, (b) Zoom-in of the anticrossing gap in (a) revealing the line [e] running 
through the anticrossing. (c) Voltage sweep along line "B" in figured! 



of a selection rule (instead of incorrectly assigning them to an additional state of the 
molecule). The intensity of these spin-forbidden transitions is proportional to |D| 4 . 

However, as the field is increased further the spin-allowed peak approaches the 
ground-state line (boundary of the SET regime) and the transport spectrum displays a 
pronounced anticrossing. This maps out the anticrossing in the N = 2 energy spectrum, 
as discussed in the introduction (cf. (OQ)). The gap between the oII/oIVb peaks directly 
gives access to the magnitude of the Dzyaloshinskii-Moriya vector: 

A\/ gap = 4=|D| (5) 
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However, the zoom in figure 6(b) reveals an additional resonance running straight 



through the anticrossing gap. This deviation from the energy spectrum is a true non- 
equilibrium transport effect. This is confirmed by our calculations including relaxation 
rates: this transition (marked by [e] in figure [5]) becomes suppressed since it starts in an 
excited state. In contrast, the anticrossing is robust against relaxation, since it involves 
transitions (marked by [d] and [f]) starting from the N — 1 ground state. 

A similar signature of the Dzyaloshinskii-Moriya vector is found, if one measures 
dl /oIVb as function of Vb and B starting from the N = 2 Coulomb blockade (CB) 
regime (line "B" in figure H|). The calculated result is shown in figure 6(c) The 



boundary of the SET regime shows an adiabatic crossover, rather than a sharp transition, 
around B z w J. In addition an anticrossing of transport excitations appears with 
conductance amplitudes which show an unusual, pronounced bias and field dependence 
due to the singlet-triplet mixing. For low field the excitation branch [g] is dominant 
since the transition is spin-allowed. As one goes through the anticrossing this branch 
becomes suppressed since it has exchanged its physical character with branch [h] which 
is spin-forbidden at low fields. We note that here the DM interaction gives rise to 
two effects: the anomalous g-factor of the spin-forbidden section of the branches and 
their anticrossing. Finally, we mention that in the presence of relaxation the branch [h] 
becomes suppressed since the initial state is an excited state, in contrast to branch [g] 
which is unaffected. 



We now discuss the effect of the magnetic field orientation presented in figures 7(a 



and 7(b) As explained in the introduction, for B || D the anticrossing in the energy 
spectrum vanishes. Therefore if the external magnetic field is rotated, the size of the 



gap in the transport spectrum will oscillate as 7(a) and 7(b) show. Interestingly, the 



non-equilibrium peak in the anticrossing gap (process [e]), shows also a clear dependence 
on the direction of the magnetic field due to the admixture between \S) and |T ) for 
B || D. However, when rotating the magnetic field, the field magnitude for which singlet 
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Figure 7. Magneto-conductance maps for regime I: dependence of dl/dVB on Vb and 
the relative angle <p between B and D for magnetic field magnitude (a) B z = 1.134 J 
and (b) B z = 1.08J. 
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and triplet should (anti-)cross is shifted. This is due to an interplay between orbital 
Zeeman effect and spin-orbit coupling leading to an anisotropic renormalization of the 
effective g-tensor in (CQ). This effect is captured by our exact diagonalization of the 
model H m . Furthermore, the orbital Zeeman effect induces oscillations in the energy 
and transport spectrum when rotating the field. Therefore the observed oscillations with 
period tt in figures 7(a)|7(b) cannot simply be attributed to the DM interaction alone. 
This underlines the importance of transport calculations using the detailed model H m . 

We now explain how the above transport signatures allow one to determine the 
axis of the DM vector experimentally. Starting with an arbitrary orientation of the 
magnetic field, one first measures the dependence on the field magnitude in search for 
an anticrossing resulting in a map similar to figure 6(a) Having found the anticrossing, 
one then fixes the field strength and rotates the field by 2ir. One then observes an 



oscillation of the anisotropic magneto-conductance similar to that shown in figure 7(b) 



with a remnant gap at the crossing points. Next one can reduce this remnant gap 
to zero by a combination of adjusting the magnetic field (to compensate the g-factor 
renormalization) and changing the plane of rotation until it includes the D-vector. 
One has determined the axis of the D vector when the measured map corresponds 
to figure 7(a) where for B || D a crossing is observed. We emphasize that the maximal 
gap observed in this figure for ip = it/ 2 does not allow the DM vector magnitude to 
be read off. Instead, one needs to orient the magnetic field perpendicular to the D- 



vector and measure the map corresponding to figure 6(a) The measured anticrossing 



gap directly gives the magnitude of the D-vector via ((HJ). We note that this last step is 
necessary due the renormalization of the g-factor. 



3.2. Ferromagnetic electrodes 

Although with the measuring scheme proposed above one can determine the axis of the 
D-vector, it is impossible to infer whether it is aligned antiparallel or parallel to this axis. 
This is due to the fact that the effective Hamiltonian (0Q) is invariant under simultaneous 
reversal of D and interchange of the spins Si and S 2 (leaving the electrodes unaltered). 
Clearly, with ferromagnetic leads this symmetry is broken and one can therefore probe 
the chirality of the ground state. Although nanoscale spin-valve structures can be 
realized nowadays (251 [2S1 [27], we note that to perform measurements with polarized 
electrodes in a magnetic field in the tesla range one may make use of alloys with large 
coercive fields and the effect of micromagnetic structure of the tunnel contact region [28J. 

To clearly probe the chirality both polarization axes of the electrodes needs to be 
perpendicular to the plane spanned by the external magnetic field and the D-vector. 
To obtain a maximal effect the field needs to be perpendicular to the axis of the DM 
vector, i.e., B _L D. Note that these requirements can be met in principle since the 
directions can be determined by first applying the measurement scheme of the previous 
section to the case ferromagnetic electrodes. We make use of the fact that the tunnelling 
onto and off the monomers is sensitive to the relative orientation of the monomer spin 
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(a) (b) 

Figure 8. Probing the chirality of the ground state with ferromagnetic leads with 
polarization p (right) and —p (left), (a) p > 0: for each monomer the projection of its 
spin onto the polarization vector of the adjacent electrode is positive and the current 
is not inhibited, (b) p < 0: reversal of the polarizations suppresses the current relative 
to situation in (a). 



to the polarization vector of the connected electrode. Therefore, for antiparallel lead 
polarizations the monomer spins in the N = 2 ground state will have orientations 
which are either both favourable (figure 8(a)) or both unfavourable for the transport 
(figure 8(b) ). The transport measured for two opposite antiparallel polarizations is thus 
expected to probe the chirality of the spin-canting due to the DM interaction. 

In figure [9] we present dl/dVs vs Vb curves calculated along line "A" in figure H] 
at the anticrossing (|B| = J) for different degrees and sign of the polarizations while 
maintaining the antiparallel configuration. The degree of polarization for the right 
electrode is defined as p = (u^ — v^)j{y^ + where a =f , j denotes the spin quantized 
in the direction of the polarization vector n of the right lead and v a denotes the density 
of states for spin a. For the antiparallel configuration considered here the degree of 




Figure 9. Conductance dl/clVs as function of Vb along line "A" in figured] for fixed 
magnetic field BID and B z sa J. The green vertical lines indicate the transition 
energies. The small shifts of the [d] and [f] peaks is due to the combined effect of 
the large difference in the tunnel rates for these states and the finite temperature (see 
[221 [S]). Note that the non-equilibrium mid-gap peak [e] discussed in section ETT1 is 
sensitive to the polarization as well. 
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polarization of the left electrode equals —p. Figure M shows that the ground to ground 
conductance peak [d] is clearly enhanced for p > and suppressed forp < 0, whereas the 
peak [f] involving the N = 2 excited state with opposite chirality shows the opposite 
behaviour. This effect is a clear signature of the definite spin-chirality induced by 
the DM interaction. The effect is most pronounced near the anticrossing where the 
spin-multiplets strongly mix due to DM interaction, and does not even require nearly 
complete polarization of the electrodes (\p\ ~ 1). We mention that even at fields below 
the anticrossing the effect can be significant if the polarization is sufficient. 

From a measurement resembling figure [9] one can experimentally conclude that the 



ground state spin structure is as shown in figure 8(a), Given the polarization n of 



the right electrode, one infers that the D vector points in the direction ±B x n if an 
enhancement /suppression of the ground to ground state peak is measured. From this 
one then finds the microscopic spin-canting Si x S2 in the N = 2 ground state, which 
is antiparallel to D in order to minimize the DM interaction energy. 



4. Conclusion 



We have predicted non-linear transport signatures of the Dzyaloshinskii-Moriya 
interaction for the generic case of a molecular dimer. We have shown how both the 
magnitude and the orientation of the vector characterizing this interaction can be 
determined from three-terminal transport measurements. The violation of the spin 
selection rules results in weak spin-forbidden electronic transitions with an apparent 
tripled g-factor in the transport spectrum. Most notably, the adiabatic change in 
the magnetization induced by the DM interaction causes a strong dependence of 
conductance peak amplitudes on the field strength and orientation. We show that when 
extracting the axis of the DM-vector from the magneto-conductance map, one has to 
account for the orbital Zeeman effect and the resulting anisotropic renormalization of 
the ^-factor. The proposed measurement scheme is robust against significant relaxation 
of the magnetic excitations. When probing the dimer with polarized electrodes with 
antiparallel magnetizations, the canted spin structure induced by the DM interaction 
leads to a strong sensitivity to the reversal of the polarization vectors. This effect allows 
the absolute orientation of the DM vector to be inferred. 

On a more general level, this work shows how transport can probe complex 
intramolecular excitations which arise when monomers hybridize and their local 
excitations become coupled. In analogy to this work, mechanical excitations of two 
coupled monomers were shown [21] to result in SET transport effects signalling pseudo 
Jahn- Teller dynamics and the breakdown of the Born-Oppenheimer approximation. For 
the three-electron dimer states an additional interplay of these vibrations with the 
monomer spins resulting from double-exchange allows for vibrational spin-spectroscopy 
and even results in a vibration- induced spin-blockade [22]. All this goes to show that 
the quantized degrees of freedom of molecular "double quantum dots" and their various 
interactions lead to a broad range of non-linear transport effects which provide detailed 
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insight and may prove useful in nanoelectronic devices. 
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